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We study the interplay between a XY anisotropy 7, exchange modulations and an external mag- 
netic field along the z direction in the XY Z chain using bosonization and Lanczos diagonalization 
techniques. We find an Ising critical line in the space of couplings which occur due to competing 
relevant perturbations which are present. More general situations are also discussed. 
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Introduction. The competition between different relevant perturbations which can render a system critical in a 
certain domain of the couplings space has been studied by many authors [1-8]. Within a low energy description of 
many different ID lattice models, one often finds the double-frequency sine-Gordon model [3], that is, a U(l) scalar 
field with two perturbations of different frequencies. In the case where both perturbations are relevant and for a 
certain frequency ratio, it has been conjectured that an Ising criticality can arise [3]. Such a situation has been found 
in the study of different lattice systems (see e.g. [6] and references therein). More recently, the so-called self-dual 
sine-Gordon model has been studied in [6]. 

A simple realization of this situation is found in an exactly solvable case, i.e. in the dimerized XY chain [8,9] and 
a qualitative analysis of the phase diagram has been presented for more general cases in a field [8]. The effect of an 
XY anisotropy has also been studied for the Fibonacci XY chain in a magnetic field where it was shown that the 
rather involved staircase structure of the magnetization curve gradually disappears by increasing the anisotropy of 
the spin exchange interactions [10]. 

In the present paper we present a unified picture of the above mentioned effects and describe, by using bosonization 
as well as numerical techniques, generic situations in which a spin gap opening mechanism [11] competes with an XY 
anisotropy to render Ising criticality. This issue was studied first in [7,8]. In [8] the interplay between a gap-opening 
perturbation and XY anisotropy was analyzed and a qualitative phase diagram was proposed. In this paper we focus 
on the Hamiltonian 

H = £ ((1 + !)SlSl +l + (1 - 7)^+1 + {S x n S* +1 + SlS y n+1 ) + AS z n S z n+1 ) -h^S^, (1) 
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which captures most of the essential aspects referred to above, and provide quantitative results which support the 
statements made in [8]. A more general situation as that arising in the Fibonacci XY Z chain in a magnetic field is 
discussed in the conclusions. 

Free Fermion Results. In preparation for the analysis of more general situations, first we discuss the XY or A = 
case which already bears some generic features. In this case the model is exactly solvable as it reduces to a bilinear 
fermionic form. As is well known, after the usual (1-d) Jordan- Wigner transformation the Hamiltonian can be readily 
diagonalized by means of a Bogoliubov transformation, leading to 
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H = ^2{E ki0 (S^h) (4, d*,o + 1 ) + E ktl (5,i,h) (4,i d M + 1 )} ( 2 ) 
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where 

£,(0,1) = "7| t 1 + cos ( fc ) + (* 2 + ^X 1 " cos ( fc )) + 2h " ± /(^> 7> k, h)] 1/2 (3) 



and 



f(5, 7, k,h) = 2 [£ 2 7 2 (1 - cos(k)) 2 + 2h 2 (l + cos(k) + S 2 (l - cos(fc)))] 1/2 . (4) 
The mean value of the magnetization is thus given by 

/7ivr\ 9{H) GS lJ^\dE k Q dE kl 1 



dh L 

k=l 



dh dh 



1 



where GS denotes the ground state, which is the vacuum of Bogoliubov particles' Fock space. (It should be noticed 
that the external field does not act as a chemical potential for these particles). 

From Eqs. (3) and (4) one easily obtains the critical field values where one of the modes becomes gapless, i.e where 
the Ising transition occurs [8] 

h c = ±v^72 } (6) 

which shows that the Ising transition can only occur for S > 7. 

The exact magnetization curve shown in Fig. 1 displays all the features we want to discuss (we focus in the region 
of h > 0). First of all, one notices that there is no actual plateau since the magnetization starts to increase as soon 
as h is turned on, and this is due to the breaking of the U(l) symmetry for 7 7^ . In the region of fields below h c 
[Eq.(6)], we are in what we call the "pseudo plateau region", where the slope of the magnetization curve is small 
and we have dominant density wave correlations. At the critical field one observes the Ising transition where the 
magnetization should behave as M — M c cx (h — h c ) (log \ h — h c \ — 1) . For h > h c , the system is in the XY phase 
where the magnetization increases more rapidly. Finally, a second Ising transition is observed before saturation, which 
occurs at h = 1 independently of the values of S and 7. For S < 7 the first transition does not occur. 

Bosonization Analysis. We discuss the results obtained so far within the bosonization approach given below, where 
we argue that the same picture is valid for arbitrary A, provided S and 7 are suitable renormalized. For small 7 and 
S one can study the effects of such perturbations using bosonization, which allows one to include A and h exactly 
through the Bethe Ansatz solution for the Luttinger parameter K. 

The large scale behaviour of the XX Z chain can be described by a U(l) free boson theory with Hamiltonian 

H = ±Jdx (vK(dJ) 2 + ^(d x <P?) , (7) 

which corresponds to the Tomonaga-Luttinger Hamiltonian. The field (j) and its dual (j) are given by the sum and 
difference of the light-cone components, respectively. The constant K governs the conformal dimensions of the bosonic 
vertex operators and can be obtained exactly from the Bethe Ansatz solution of the XXZ chain (see e.g. [12] for a 
detailed summary and references therein). We have K = 1 for the SU(2) symmetric case (A = 1) and is related to 
the radius R of [12] by K~ l = 2ttR 2 . In (7) v corresponds to the Fermi velocity of the fundamental excitations of the 
system. In terms of these fields, the spin operators read 



1 d x <j> + a : cos(2fe F x + A/27T0) : +-^- , (8) 



V2tt 

S± ~ (-1)* : e ±i ^ (bcos(2k F x + V^) + cj : , (9) 

where the colons denote normal ordering with respect to the ground state with magnetization (M). The Fermi 
momentum kp is related to the magnetization of the chain as kp — (1 — (M))tt/2. The effect of an XYZ anisotropy 
and/or the external magnetic field is then to modify the scaling dimensions of the physical fields through K and 
the commensurability properties of the spin operators, as can be seen from (8), (9). The constants a, b and c were 
numerically computed in the case of zero magnetic field [13] (see also [14]). 
The bosonized Hamiltonian including the perturbations then reads 

Hij OS = H + Ai J dx COs(V27T(j)) + A2 J dx cos(VStt4>) (10) 

where Ai ex (5, A2 oc 7. 

The scaling dimensions of the perturbations in (10) are K/2 and 2/K respectively, which in the XX case (K = 2) 
are both equal to unity. This allows one to understand within this approach the appearance of the Ising transition [6] 
since in this self-dual case and for h = one can map the bosonic system into two Major ana fields, and the critical 
line is given by Ai = ±A2 where one of the masses of the two Majorana fields vanishes rendering criticality. As soon 
as the external field is turned on the masses of these Majorana fermions change and we then need a bigger value of 
Ai (i.e. S) to find the transition [8] [see (6)]. 

This can be extended for A > 0, since its only effect is to modify the scaling dimensions of the two perturbations 
through K(A,h). It is interesting to note, however, that for A > the cosine of the dual field is less relevant, 
becoming barely marginal for A = 1 and h = (if (1, 0) = 1), opening the question if one should still expect the Ising 
transition to occur. One can argue that this is so by analyzing the one-loop renormalization group (RG) equations, 
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where it can be seen that due to the renormalization of K, which (up to this order in the RG) flows to the XX value 
at large scales for certain initial values of the couplings, the competition between the two cosines can still happen 
giving rise to the Ising criticality even in this limiting situation. From this qualitative analysis one expects a critical 
line for h = given by 

8(xY (11) 

wit h y = K * f l£~*'tf ,A) and K eff is the large scale value of K. As mentioned above, for < A < 1 , the one loop 
RG equations indicate that K e f f = 2 , therefore implying that v = 1 . A thorough RG analysis should be performed 
nevertheless, since this large scale value of K e ff is attained in a strong coupling regime, where the analysis ceases to 
be valid, but this is out of the scope of the present paper. We provide numerical evidence for the validity of eq. (11) 
with v = 1 below. 



Numerical Results. In what follows, we examine this latter conjecture by studying numerically the extreme SU (2) 
case (A = 1). On the other hand, this enables an independent test of the bosonization scenario within non- perturb at ive 
regimes. 

First, we computed the gap of Hamiltonian (1) at (M) = by means of Lanczos diagonalizations [15] of finite 
chain lengths L with periodic boundary conditions. If L/4 is an integer it turns out that the ground state is even, i.e. 
exp (in ct+<t~) = 1 , otherwise the gap spectrum has to be computed within the odd subspace of H [16]. In Fig. 
2 we display the results so obtained for 14 < L < 22 (see [17] ). As expected, finite size effects become noticeable on 
approaching the critical regime encompassed within the pronounced gap minima, whereas convergence towards the 
thermodynamic limit becomes typically logarithmic [18]. Interestingly, the locations of these minima result however 
fairly independent of the system size, thus facilitating the estimation of the critical line. For L < 12 however, the 
spectrum gap increases monotonically with |7|. Therefore, to avoid misleading results we dismiss the data of those 
small sizes. In turn, this impedes the usage of most of the recursive extrapolation algorithms encountered in the 
literature [18], which show their full strength only if the can be iterated several times on a large sequence of sizes. 
Thus, we content ourselves with a standard (logarithmic) gap extrapolation of the form g ~ gap (L) + A/L B . Though 
semi-quantitative [i.e. A = A(^y, S) , B = £(7, 8) ], the results shown in the inset of Fig. 2 nevertheless indicate a gap 
closing at the minima of the finite size data. 

In studying numerically the massless regime obtained from such minima however, it should be borne in mind that 
finite size corrections to the gap of the homogeneous Heisenberg chain (7 = £ = 0), vary slowly as ln(lnL)/ In 2 L , 
thus affecting the results over a wide range of sizes [19]. Therefore, we are confronted to restricting considerations to 
regimes where either I7I or \S\ are not too small. Fig. 3 exhibits our estimation of the critical line down to S = 0.15 
above which however, a wide linear regime shows up, in agreement with the bosonization approach. Similar results 
were observed for < A < 1 as well. 

Secondly, we turn to the GS magnetization curves of Hamiltonian (1). These were calculated numerically as 

^2j(GS\ (j z - \GS) , where in general the parity exp (in ^Z n o~^cr~) °f \GS) comes out to be both field and size 
dependent. In spite of the existence of a whole massive regime away the critical line, the magnetization grows linearly 
for small applied fields h (a feature which holds also for A = 0), rather than displaying a plateau at (M) = 0. 
This is shown by the inset of Fig. 4(a). Since this latter regime is dominated by a finite gap width (i.e. short range 
correlations), finite size effects become negligible there. This renders our low field results most reliable, until the 
first parity change in \GS) occurs (signaled by an abrupt increase of the magnetization), probably associated to 
the emergence of gapless modes such as those referred to in the free fermion case. When selecting the Hamiltonian 
parameters in a massless (zero field) point, the magnetization still grows linearly with h though with a substantially 
bigger amplitude, as is shown by the inset of Fig. 4(b). This results in the removal of all the pseudo-plateaux observed 
in Fig. 4(a) with a massive 7 — S point. Also, notice that near the brink of saturation in either case the susceptibility 
tends to diverge, alike the XY situation. 

Finally, we address to the predicted Ising behavior of the magnetization curves near the critical fields h c discussed 
in the previous sections. Upon estimating the former with the fields h c (L) involving the first parity jump in the 
GS of a finite chain, we obtained a fair logarithmic regime actually applying over more than two decades in (h — 
h c ) (log \ h — h c \ — 1) . This is displayed in Fig. 5. The value of the upper line slope results almost independent of the 
system size and is indicative of the logarithmic regimes entailed both by the free fermion and bosonization approaches. 

Conclusions. To summarize, we have studied how the tendencies towards the formation of massive spin excitations 
(through dimerized ^-exchanges), and towards the XF-ordering (via pairing of 7-interactions), compete with each 
other. In the bosonization picture, these tendencies manifest themselves in the existence of two (competing) relevant 
interactions and bring about the Ising transition, which in fact is connected to the same one driven by an external 
field. Due to the breaking of the U(l) symmetry, the latter does not couple to any conserved quantity and therefore 
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the magnetization process gets essentially modified with respect to the 7 = case, i.e. it ever increases in spite of 
the presence of massive regimes. 

Both the bosonization and the numerical analyses support the picture that the physical mechanism rendering the 
Ising transition is quite general and valid both in the strong and weak coupling limits. Following the ideas given in 
[20] one can conjecture that a similar picture for each plateaux in the XYZ Fibonacci chain holds. The mechanism 
in that case is the same, where the relevant operator coming from dimerization is replaced for each plateaux by the 
operator commensurate at the corresponding frequency of the Fibonacci Fourier spectrum. From this one concludes 
that a similar picture as that found in the XY case [10] is also valid for generic XYZ Fibonacci chains. 

The numerical analysis lent further support to the bosonization results within several non-perturbative situations. 
We found evidence of a massless line along with the expected Ising like behaviour near critical fields. In turn, this 
suggests that most basic features of the fully interacting system can be captured by the free fermion picture discussed 
above. However, the issue as to whether or not all of the GS parity changes induced by the field [Fig. 4(a) ] correspond 
to Ising transitions in the thermodynamic limit, remains quite open. 
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FIG. 1. Exact magnetization curve of the XY case. The 
inset displays the log arithmic singularity dominating the sus- 
ceptibility at h c = \J S 2 — 7 2 . 
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FIG. 2. Gap spectrum for different sizes of Hamiltonian 
(1) with A = 1 . The inset exhibits the (estimative) gap 
extrapolations to the thermodynamic limit. 
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FIG. 3. Gap less line of Hamiltonian (1) in the (j,S) cou- 
pling parameter space arising from the positions of the gap 
minima of finite samples with A = 1 and (M) — . 



FIG. 4. Magnetization curves for A = 1 . (a) The 
piece-wise continuous behavior is related to parity changes 
in the ground state, (b) Results arising for a massless (7,^) 
point. In each case, the inset exhibits the linear response at 
low- field regimes. 
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FIG. 5. Magnetization behavior for A = 1 near critical 
fields referred to in the text. The slope value suggests the 
emergence of a logarithmic regime as that observed in Fig. 1. 
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